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Abstract 

We study the three-quark light-cone amplitudes of the proton including quarks' transverse mo- 

', menta. We classify these amplitudes using a newly-developed method in which light-cone wave 

functions are constructed from a class of light-cone matrix elements. We derive the constraints 

on the amplitudes from parity and time-reversal symmetries. We use the amplitudes to calculate 

CN ■ the physical observables which vanish when the quark orbital angular momentum is absent. These 

O 
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include transverse-momentum dependent parton distributions Aqx(x,k±), qT( x ,k±), 5q(x,k±), 
and 6qL(x,k±), twist-three parton distributions gT{x) and hi(x), helicity-flip generalized parton 
distributions E(x,£ = 0,Q 2 ) and its associates, and the Pauli form factor F2(Q 2 ). 
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I. INTRODUCTION 



In the simplest model of the proton, the valence quarks are moving in the s-orbit and the 
proton spin is constructed from the quark spins. This picture has been ruled out definitely 
by the EMC data Q and the follow-up experiments, SMC, E142, E143, E152, E153, and 
HERMES 0. One then expects that the quarks in the proton must have non-zero orbital 
angular momentum and/or the gluons must carry part of the proton spin. When the orbital 
angular momenta of quarks and gluons are none-zero, the nucleon is intrinsically deformed. 
An interesting question is how big is the deformation? or how much does the quark orbital 
angular momentum contribute to the spin of the nucleon? A few years ago one of us derived 
a sum rule, relating the total angular momentum carried by quarks to an integral over the 
relevant generalized parton distributions (GPDs)|3|, [|]. The orbital angular momentum can 
be obtained by subtracting the measured quark helicity contribution. This observation has 
lead to very active studies of GPDs and their measurements in hard exclusive processes^], 

B0. 

There are, in fact, many observables which are potentially sensitive to the orbital an- 
gular momentum of the quarks, although they do not directly measure the orbital angular 
momentum itself. For example, the Pauli form factor -^(Q 2 ) of the proton, which has 
been measured recently using the recoil polarization technique at JLab || ||, has long 
been recognized as an observable related to the light-cone amplitudes with non-zero or- 
bital angular momentum [10[ [XT , [12]. The twist-three spin-dependent structure function 
gr{x) = gi{x) + g2{x) vanishes in the naive parton model, and its interpretation requires 
quarks' orbital motion [fl3|, |i~4| , The generalized parton distributions such as E(x,l;,Q 2 ) 
were introduced to characterize quark orbital angular momentum in the first place |3], [16], [17| . 
And finally, the transverse-momentum dependent parton distributions 
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contain rich and direct information about the quark orbital motion [24, 25, 26, 27]. A nat 



ural question is can one correlate all these observables in terms of light-cone wave functions 
in a minimally-model-dependent way? This paper is an attempt in this direction. 

Our approach starts with the light-cone expansion of the proton state truncated to the 
minimum Fock sector with just three valence quarks |2S, [ZIJ. Unlike the light-cone ampli- 
tudes studied in the exclusive processes J5D], |3T], [J5[ [33|, 0, we keep the full dependence on 
the transverse momenta of the quarks. Some previous studies in this direction can found in 
Refs.|35, 36, 37|. Clearly, this truncation is ideal: we do not have the gluon and sea quark 



degrees of freedom; we do not have manifest rotational symmetry and gauge invariance. 
However, the truncation can be improved systematically by including higher Fock compo- 
nents, for example, adding one or more pair of sea quarks and/or one or more gluons, etc. 
Therefore, in this study we will ignore certain artifacts arising from this specific truncation. 

Based on our work here, one can go on to parameterize the three-quark amplitudes and 
fit to the experimental data. The amplitudes thus determined are phenomenological, and in 
principle cannot be compared with those solved from QCD directly. However, a comparison 
helps to determine the importance of the three quark amplitudes vs. multi-particle Fock 
amplitudes. One may also consider the phenomenological amplitudes as the effective proton 
wave function after integrating out the gluon and sea quark degrees of freedom, although 
the integrations are hard to implement in practice. 

By committing ourselves to the light-cone amplitudes, we are also committing to the 
light-cone gauges A + = [38]. O ne can in principle add light-cone links between the quark 
fields to make the light-cone amplitudes gauge- invariant, but we do not find a compelling 
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reason to do this. One subtlety about the light-cone gauge is that it requires additional 
gauge fixing [39, 50, fl]]. Depending on whether the additional gauge condition is time- 



reversal invariant or not, the wave function amplitudes are real or fully complex. In the 
latter case, the final state interaction effects are contained in the amplitudes [[27], |33 



Our plan of the presentation is as follows. We start in Sec. II by writing down the matrix 
elements of three-quark light-cone operators between proton state and the QCD vacuum, 
which serve to define a complete set of light-cone amplitudes within the truncation. We 
simplify these matrix elements using color, flavor, spin, and discrete symmetries. We then 
invert them in Sec. Ill to find the three-quark light-cone wave functions. As we have 
explained in Ref . [ |4|] , there are many advantages in this construction. One is that the cut- 
off dependence of the light-cone amplitudes can be studied directly from the matrix elements. 
We derived the scale evolution equation of the individual Fock contributions to the parton 
densities in [ 44[ . Another is that these matrix elements may be calculated in the nucleon 



models or lattice QCD: it is not necessary to construct a model of the nucleon directly in the 
light-cone frame. Finally, the large-momentum behavior of the nucleon amplitudes can be 
determined from Bethe-Salpeter-like equations governing the matrix elements. In Sec. IV, 
we calculate a number of physical observables which are directly linked to the amplitudes 
with non-zero quark orbital angular momentum. 



II. THREE-QUARK LIGHT-CONE MATRIX ELEMENTS 

In this section, we are interested in classifying the matrix elements of three-quark oper- 
ators between the proton and QCD vacuum. Since much of the discussion depends on the 
use of light-cone coordinates and light-cone gauge, we start by reminding the reader some 
salient features of the light-cone technology, and a more detailed exposition can be found in 
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The light-cone time x + and coordinate x~ are defined as x ± = l/^/2(x° ± a; 3 ). Likewise 
we define Dirac matrices ^ = l/v^2(7° ±7 3 )- The projection operators for Dirac fields are 
defined as P± = (1/2)7 =F 7 ± . Any Dirac field ip can be decomposed into ip = ip + + ^_ with 
ip± = P±ip. ip+ is a dynamical degrees of freedom and has the canonical expansion, 



-i(fc+r-fe_L-H) 



+d[(k)v(k\)e i{k+ ^-^ ) ] . (1) 

Likewise, for the gluon fields in the light-cone gauge A + = 0, A± is dynamical and has the 
expansion, 

A 

+a{(k)e*(k\)e i ^~- U ^ ± A . (2) 
and A~ are dependent variables. 



The key observation in Ref. is that the light-cone Fock expansion of a hadron state 
is completely defined by the matrix elements of a special class equal-light-cone-time quark- 
gluon operators between the QCD vacuum and the hadron. These operators are specified 
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as follows: Take the + component of the Dirac field and the + _L-component of the 
gauge field F +± . [We sometimes label the _L components with index % = 1,2.] Assume 
all these fields are at light-cone time x + = 0, but otherwise with arbitrary dependence 
on other spacetime coordinates. Products of these fields with the right quantum numbers 
(spin, flavor, and color) define a set of operator basis. [This has been done in the past for 
light-cone amplitudes in which all fields are separated along the light-cone, see for example 
3tJ.] Clearly, these operators are not gauge-invariant although one can gauge- invariantize 



them by inserting light-cone links extending from the locations of the fields to infinity for 
every field. Since all fields are at different spacetime points in the transverse direction, the 
operators do not require additional renormalization subtractions apart from the usual wave 
function renormalization. The matrix elements of all these operators between the hadron 
state and the QCD vacuum contain complete information about the hadron wave function. 

In the following two subsections, we will use the above method to study the matrix 
elements of operators constructed out of three quark fields between proton and the QCD 
vacuum. We discuss separately the cases for the two up-quarks coupling to helicity-zero and 
one. 



A. Two Up-Quarks Coupling To Helicity-Zero 

When the two up-quarks are coupled to helicity-0, we define the following two matrix 
elements, 

(0|^ a (ei)C7 + <(6)^(e 3 )e afcc |P) 
= 0«(1, 2, 3) T5 t/ + + {id\^ 2 \l, 2, 3) + z^< 2 >(2, 1, 3))y B ?U+ 
+d[d^ 3 \l,2,3)U + , 

(3) 

(oi^ a (ei)c 7 + 75<(6)^(e3)6 a6c b) 

= 0^(1, 2, 3)U + + (idi^ 5 \l, 2, 3) - ^0 (5) (2, 1, 3))f U + 

+^ 75 (6) (l,2,3)f/ + , (4) 

where a,b,c = (1, 2, 3) are color indices of the quarks, C is the charge conjugation matrix 
(C = ry 2 7°), and \P) is the proton state with momentum P M and Dirac spinor U+(P). The 
arguments 1, 2, and 3 in the scalar functions <fy^ stand for (^2~,^2±), and £^~,£3_i_), 

and the functional dependence on transverse coordinates are of type ■ only. The up 
and down quark fields are represented by u and d, respectively. The index i goes over the 
transverse coordinates 1, 2. In the above equations, we have used the symmetry relations 
(T stands for transpose) 

< a (£i)<?7 + <(£2)e a6c = ^ a (£ 2 )C 7 + <(6K fec , (5) 



< a (£i)C 7 + 75<(6K bc = -< a (£2)<?7 + 75<(£i) 



e 



which imply 



0«(1,2,3) = 0«(2,1,3), 0( 6 )(1,2,3) = ^(2,1,3) , 

0^(1, 2, 3) = -0 (3) (2,1,3) , (4) (1,2,3) = -0^(2, 1,3). (6) 
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The matrix elements have an overall dependence on £1 + £2 + £3 through a phase factor 
e H6H-&+&)-- p / 3 ) which we assumed has been factorized implicitly. The remaining dependence 
is on the differences of the coordinates, for example, £1 — £2 and £2 — £3, because of the 
translational invariance. Therefore, the partial derivative d\ is not independent of d\ and 
&2, and has been omitted. We have not included the structure d l Y in these equations because 
it is the same as 9*7*75 when acting on U+ using = 0. We use the shorthand 

qi _ e tjgj ^12 _ Similarly the structure d[dld{'j : 'U + can be reduced to the existing ones. 

We can simplify the above equations by adding them and applying the chiral projection 
operators P L>R = (1 =p 7s)/2 

2(0|^(6)C7 + < i (e2)^(6)6 o6c |P) 
= 2, 3) - 0^(1, 2, 3)] U +R + [0^(1, 2, 3) - 0^(1, 2, 3)] U +R 

+ [id\ (0^(1, 2, 3) - 0^(1, 2, 3)) + id\ (0^(2, 1, 3)) + 0^(2, 1, 3))] 7^+L • (7) 

The two-up quarks have been paired to helicity-zero, and the remaining down quark is 
projected to be right-handed. We can further simplify the above equation by introducing 
the new amplitudes, 

(0|^(6)C7V +L (6)^(6)e afe 7v / 6|P) 
= 2, 3)U +R + i 2 d{idi^ 2 \l, 2, 3)U +R + [id\^\l, 2, 3) + id^(l, 2, 3)] YU +L (,8) 

where -0's have no specific symmetry properties because the two up quarks are not in the 
same helicity state. If the nucleon is right-handed, we split the above equation into two 

abc r -, 

(0\ul R (^)C 7 +ul L (^)dl R (Cs)-^\P "{) = ^«(l,2,3) + z 2 aiz^ 2 )(l,2,3)J U +R , (9) 

,abc 



In our convention d % = —d/dx l and 7* is that of Bjorken and Drell [H5 



B. Two Up-Quarks Coupling To Helicity-One 



When the two up quarks are coupled to helicity- 1, we define the matrix element, 
(0\uf(^)Cta +t u^ 2 )d c + (^)e abc /VQ\P} 
= ^(1, 2, 3)7*75 14+ (^i0( 8 )(l,2,3) + ^^ 8 )(2,l,3)) 7 5f/+ 
+ (^(1,2,3) + «9 2 : (9) (2,1,3)) U + 

"zc^V (10) (i, 2, 3) + ^z^V (10) (2, 1, 3) 



+ 



+id[ i idiU {11) (1, 2, 3) + i 2 {d[di - c^')0 (12) (1, 2, 3) 



7 J 75^- 



+ 5 



(11) 



where the parentheses on a pair of indices indicate symmetrization and subtraction of the 
trace. The symmetry between two up-quarks yields 



(7) (1,2,3) = (7) (2,1,3) , 0( n )(l,2,3) = 0( n )(2,l,3) 
0( 12 >(1,2,3) = -0( 12 )(2,1,3). 



(12) 
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When inserting an additional 75 between the two up quark fields, and using cr +tr y 5 
— ie u cr +jf , we obtain a similar expression, 

(0K(6)Cm +i 7 5<(6)^(6)e a6 7V6|P) 

= 0^(1, 2, 3)7*17++ (^(1,2,3) + ^ 8 )(2, 1,3)) 75 tf + 

+ (^W(1.2,3)+^ (9) (2,1,3))?7 + 

ili^'V 10 ^!, 2, 3) + z,9^V (10) (2, 1, 3) 



— 1 



+ <^ 2 V (11) (1, 2, 3) + i\d\di - 8^)^(1, 2, 3) ^U. 



'+ 1 



(13) 



where we have used 7* = —ie^rf^s, valid when acting on U + . 

Combining Eqs. ( |TT|) and (|T3|) , we obtain the matrix elements of the quark fields in 
definite chirality states. For example, when the two up-quarks are right-handed, and the 
down quark left-handed, the total helicity of the quarks is 1/2. In this case, we have 

(01^(^)^+^(6)^(6)^7 v^lP T) = [0 (7) (l,2,3) + z(za 1 z4)0 (12) (l,2,3)] yu +R , 

(14) 

where we have used the relations YU+r = —ij l U + R and d^^U R+ = with d\ = 
d\ — ie^d{. The isospin symmetry imposes the following constraint, 

+ (0|<«(6)«^<K(6X i fe)e abc |P T) 

+ (0|^(6)C^+ i < ii (6)« +i (e3)e a6c |7 J T) = . 
Using Fierz identity, one can show 

ul R ^i)Cia + % R ^ 2 )u c +L (^)e abc = - M ^(^)C7 + < fi (6)7^ + L(6) e abc , 



(15) 



di R {^)Cia^u b +R ^ 2 )u\ L {^)e 



abc 



< a R (6)C , 7 + < iJ (e 3 )7^ + L(6) 



.abc 



(16) 



which means that </>( 7,12 ) are not independent amplitudes, 

(7) (1,2,3) = V (1) (l,3,2) + ^ (1) (2,3,1) , 
(12) (1,2,3) = ^(1,3,2) -^(2, 3,1) . 

Hence we write, 



(17) 



.abc 



^W +R (ii)Cia +l u\ R (i 2 )dX L (i z ) — \P t) 

3, 2) + V (1) (2, 3, 1) + • (iftifc) (^(1, 3, 2) - ^(2, 3, 1)) 



7*C7. 



+R 



(18) 



which contains no new amplitude. 

If all quarks are right-handed, the total helicity of the quarks is 3/2. We need one unit 
of orbital angular momentum (projection on the z-direction) to construct the helicity of the 
proton. We write, 



2(0|^(6)« ( x +i < i? (6)^fe)e afe 7v / 6|P T) 
[idt^ih 2, 3) + z<r</> (8+9) (2, 1, 3)] U +R , 



(19) 
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where 0( 8+9 ) = 0( 8 ) + 0( 9 ). The isospin symmetry leads to the following relation 

^ 8+9 \l, 2, 3) + ^ 8+9 \l, 3, 2) - 0( 8+9 )(3, 1,2)- 0( 8+9 )(3, 2, 1) = . (20) 

If we define 

^ 8+9 )(l,2,3)/2 = ^ (5) (l,2 ! 3)-^ (5) (l ) 3,2) , (21) 
the above constraint is solved. Therefore, 

,abc 

(0|^(6)^ + ^(6)^fe)^|P T) 

= [id[- (V (5) (l,2,3) -^ (5) (1,3,2)) +idt (^(2,1,3) -^(2,3,1))] U +R , (22) 

which defines a new amplitude tp( 5 \ 

When the two up quarks are left-handed and the down quark right-handed, the total 
quark helicity is —1/2. We again need one unit of orbital angular momentum to construct 
the proton helicity. The orbital angular momentum can either be on the first or second 
particle, 

2(0|^l(6)Cm +l < L (6)^fe) e a6 7v / 6|P T) 

= ^i + 0( 8 - 9 )(l,2,3) + ^ + (8 - 9) (2,l,3)]f/ +i? , (23) 

where d[ + = d\ + ie^d{ and so on. Using the constraint from the isospin symmetry, one can 
express 0^ 8 ~ 9 ^ in terms of and 

^(l, 2, 3)/2 = V (4) (3, 1, 2) - ^(3, 1, 2) - ^(3, 2, 1) . (24) 

Hence we write 

f abc 

(0|< a L (6)Cm + V +L (6)^(6)^|P T) 

= [id\ + ty«(3, 1,2)- ^ (3) (3, 1,2)- ^ {3) (3, 2, 1)) + (1 <-> 2)] U +R , (25) 

which does not contain any new amplitude. 

When all quarks are left-handed, the total quark helicity is —3/2; two units of orbital 
angular momentum are need to construct the nucleon helicity. Taking the appropriate 
projection of the above equations, we get 

2(01^1(6)^^^(6)^(6)6^7 V6|PT> 

id[ i+ idiU (10) (l, 2, 3) + ^ + ^ ) {1O) (2, 1, 3) + id[ i+ idiU (U) Ci-, 2, 3)1 i 3 U +R , (26) 

where ( ^ lo - ) describes the state in which the same quark carrys two units of angular mo- 
mentum, and describes each up quark carrying one unit of angular momentum. The 
constraint from the isospin symmetry relates to c^ 10 ** 

(11) (1, 2, 3) = (1O) (1, 2, 3) + (1O) (2, 1, 3) - (1O) (1, 3, 2) - 0^(2, 3, 1) . (27) 
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Following; our convention, we rename 0( lo )(l,2,3)/2 to ^( 6 )(1,2,3), and Eq. (]2§) becomes 



^abc 



(0|^ a L (6)«a +? < L (6)^Lfe)^|P T) 
= z<9f + ^V 6) (l,2,3) + z^ + z^V (6) (2,l,3) 

+zef f ze^ (V (6) (l, 2, 3) + ^(2, 1, 3) - V (6) (l, 3, 2) - ^(2, 3, 1))] 7 J f/+H , (28) 

which contains the new amplitude ^M 6 ). 

III. THE THREE-QUARK LIGHT-CONE WAVE FUNCTION OF THE PROTON 

The matrix elements considered in the previous section can be calculated in models such 
as the MIT bag or QCD sum rules or lattice field theory. Once they are known, one can 
invert the equations to get the light-cone wave function of the proton in the minimal Fock 
sector. In this section, we work out the light-cone wave function in terms of these matrix 
elements. The expression allows us to calculate the physical observables of the proton in the 
next section. 

Consider first the case in which all quarks are coupled to helicity 1/2. Define the measure 
for the quark momentum integrations, 

41M2M3] = V2 dxidX2dx * d2 ^Af^ 

M M U ^/2x l 2x 2 2x 3 (2tt) 9 

x2nd(l -x 1 -x 2 - x 3 )(27r) 2 (5 (2) (A?i± + k 2± + k 3± ) . (29) 

where Xj are the fraction of the nucleon momentum carried by the quarks, and ki± are their 
transverse momenta. We introduce a Fock component for a right-handed proton, 

|PT)i/ 2 = J d[l]d[2]d[3] (^ 1 )(1,2,3) + ^(^A;|-^)^ 2 )(1,2,3)) 

^abc , . 

x— u t T (i) (^( 2 )4 T (3) -4(2)^(3)) |0> , (30) 

where ip^ 1 ' 2 ' are functions of quark momenta with argument 1 represents x% and fcu_ and so 
on. The dependence on the transverse momenta is of form kij_ ■ kj± only. The ip^ amplitude 
has a pre-factor k\k\ — kfk% which signals the contributions from quarks with non-zero 
orbital angular momentum, although all magnetic quantum numbers sum to zero. Using 
the above state, the matrix elements shown in Eqs. ([5]) and (|I8|) can be calculated. We find 
the relation between the wave function amplitudes and the matrix elements 

V> (a) (£i>6:6) = / d[l]d[2]d[3]Vx 1 2x 2 2x 3 e-^^r-^-Hx) 

xe -i(x- 2 P+5 2 --fe 2 x-€2x) e -^3P + C3"-fc3X-r3x)^(«)( 1)2 ,3) , (31) 

where a = 1, 2. 
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Turn to the case where the quark helicities sum to —1/2. The corresponding Fock com- 
ponent can be written as 

|PT)-i/2 = J d[l]d[2]d[3] ((A£ + iA^ (3) (l,2,3) + (A£ + ^f)V^(l,2,3) 

abc , . 

x— (^(1)^(2)^(3) - d+ T (l)«t i (2)«t l (3)) |0) . (32) 

The matrix elements in Eqs. fllOl ) and (^) can be calculated, and the resulting f// 3,4 ) 
are related to the above amplitudes in the same way as in Eq. (]3Tf) . One might suspect 
if additional independent amplitudes can be constructed by adding terms with structure 
k x k\ — k\k x . A careful examination indicates that they can be reduced to the already 
existing ones. 

When the quark helicity is added to 3/2, the Fock component is 
|PT>3/2 = J d[l]d[2]d[3] (-^ + 1^)^(1,2,3) 

^abc , . 

x— «t T (i) (4 T (2)4 T (3) - <(2)4 T (3)J |0) . (33) 

Calculating the matrix element in Eq. (|22|) , we find t/>^(£i, £2, £3) is the Fourier transfor- 
mation of ijj( 5 \ki, &2, ks). 

Finally, we consider the case when the quark helicity adds to —3/2, the Fock component 

is 

|Pt)-3/ 2 = J d[l]d[2]d[3] (A? +i*?)(*J + i*if)^W(l l 2, 3) 

X 7I^ (1) K (2)u ^ (3) ~ u l^) d 4( 3 )) 1°) • ( 34 ) 

Using this, we calculate the matrix elements in Eq. (p8|) and find ^M 6 ^ (£1,^2, £3) is just the 
Fourier transformation of ^ 6 '(ki, k2, k^). 

The complete three-quark light-cone Fock expansion of the proton has the following form, 

\P T) = \P T>-3/2 + \P 0-1/2 + \P T)l/2 + \P T>3/2 • (35) 

Many interesting proton observables can be calculated using the above wave function as we 
will show in the next section. 

Using the results in the previous section, one can also construct a proton state with the 
negative helicity \P [). All the wave function amplitudes are the same, except that the 
quark helicities are flipped, k x ± ik y becomes k x =|= ik y , and some signs are added, 



\P[)-i/2 = j d[l)d[2]d[3) ^«(l,2,3)+i(Arf^-W)^ (a) (l 5 2,3)^ 

,abc . . 

x— ^(1) (^(2)^(3) - 4t( 2 )<(3)) |0> , (36) 
\P[) l/2 = J d[l]d[2]d[3] ((^-^)^ 3 )(l,2,3) + (^'-^)^ 4 )(l,2,3)^' 

^abc . . 

X 7I KW^t( 2 )4t(3) - 4 4 (1)<(2)<(3)) |0) , (37) 
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|P|)_3/2 = / d[l]d[2]d[3] (-)(^ + ^)V^(l,2,3) 

^abc , . 

X ^|4|(l) (4|(2)<(3) - <(2)<(3)) |0) , (38) 

|P j> 3/2 = J d[l]d[2]d[3] (-){kl - ik\)(kl - iklW\l^3) 

x^! T (l) (4 T (2)^ T (3) -4 T (2)4 T (3)) |0) . (39) 

The same expressions can be obtained from Jacob and Wick's method [Q . 

(— 1) S ~ A |P — A) = Y\PX) , (40) 

where Y is the parity operation followed by a 180° rotation in the ^-direction, and s is the 
spin and A the helicity. 

A. Constraints From Time-Reversal And Parity 

Consider the proton state |P |). Under parity, the 3-components of the proton momentum 
change the sign: (P°,P) becomes (P°, — P). Under time-reversal, the momentum changes 
in the same way. Therefore under combined time-reversal and parity, the momentum of 
the proton does not change. Neither do the quark momenta. On the other hand, helicity 
changes sign under the combined parity and time reversal, and hence the proton state |P f) 
becomes |P \). Similar changes occur for the quark states. 

Under time-reversal, one must replace C-numbers with their complex conjugates. Thus 
under the combined time reversal and parity, a positive-helicity proton state becomes a 
negative helicity one |P j) with wave function amplitudes complex-conjugated, quark helicity 
flipped, and k x ± ik y becoming k x =p ik y . Moreover, there is a sign change of (— 1) S ~ A from 
time-reversal. The resulting state |P j) is exactly the same as that in Eqs. (36-39), except 
all amplitudes are complex conjugated. Hence if time-reversal symmetry applies, all the 
wave function amplitudes must be real, 

= i = 1,2,3,4,5,6, (41) 

which is true in many model calculations. 

The above result, however, is only correct if the gauge condition is also invariant under 
the discrete symmetries. In light-cone gauge, A + = does not fix the gauge completely, 
additional boundary conditions for the gauge field must be specified |39|, ||| |44j . This addi- 



tional gauge fixing might not be invariant under the combined parity and time- reversal. For 
example, the advanced and retarded boundary conditions transform into each other under 
time-reversal; however the principal value prescription is self-conjugating. In the former case, 
the wave function amplitudes are no longer constrained by the above condition. Therefore, 
although the final-state-interaction gauge links vanish in the momentum-dependent parton 
distributions in the light-cone gauge, the imaginary part of the wave function amplitudes 



can reproduce the final state interaction effects J24], P7fl . In the following discussion, 

we use the advanced boundary condition and hence the imaginary parts of the amplitudes 
do not vanish. 
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IV. OBSERVABLES DEPENDING ON QUARK ORBITAL ANGULAR MOMEN- 
TUM 



As we have explained in Introduction, the quark orbital angular momentum in the proton 
is certainly non-zero. To determine it experimentally, one has to measure the relevant 
generalized parton distributions f3|. On the other hand, there are many other nucleon 
observables which are sensitive to the quark orbital angular momentum, although they do 
not measure it directly. The most familiar example is the Pauli form factor i^Q 2 ) of the 
nucleon. In light-cone quantization, i^Q 2 ) is an helicity-flip observable which depends on 
the interference between the wave function amplitudes differing in one unit of orbital angular 



momentum 1C, 11 



We have to caution, however, that when the gauge fields are taken into account explicitly, 
some orbital angular momentum effects discussed here cannot be cleanly separated from the 
gluon contributions. In fact, the quark angular momentum 

L q = Sx^x x (-iD)i/> , (42) 



does contain the gluon contribution through the gauge potential in the covariant derivative. 
In the example of the electron magnetic moment in quantum electrodynamics, the photon 
Fock component plays an important role flnj. Hence, some results in this section can only 



make sense if considered as an effective description after integrating out the gluon degrees 
of freedom. 

In the quark models, the need for a quark orbital angular momentum can sometimes be 
avoided by large constituent quark masses. The quarks can have significant spin flips through 
the mass term. Then some observables discussed in this section are directly proportional 



to the constituent masses In QCD, because the light flavors have negligible masses, 

the quark spin flips cannot occur and the orbital angular momentum and gluon angular 
momentum are thereby essential. 

We start by calculating the hadron helicity-flip transverse-momentum dependent parton 
distributions. The twist-three distributions are found to be simply the integrals over the 
appropriate transverse-momentum dependent distributions. After this, we calculate the gen- 
eralized parton distribution and electromagnetic form factors, with E(x,l;,Q 2 ) and F 2 (Q 2 ) 
as examples. 



A. Transverse-Momentum Dependent Parton Distributions 



The transverse momentum dependent parton distributions were first introduced by 
Collins and Soper to describe Drell-Yan production fl8|, |H and later by Sivers to describe 

A classification of the lead- 



the single spin asymmetry in hadron-hadron scattering |19], 



22 



Recently, there 



ing distributions in terms of spin and chirality can be found in Ref. 
has been much interest in measuring these distributions in the semi-inclusive deep-inelastic 
scattering ||, g|, H, (|| 0, 0, 

Let us first recapitulate the classification, and in the process, we introduce some new 
notations which we believe are easy to systematize. For an unpolarized nucleon target, one 
can introduce the unpolarized quark distribution fi(x, k±) = q(x, k±) and time-reversal odd 
transversely-polarized quark distribution h^(x,k±) = Sq(x,k±) arising from the final-state 
or initial-state interactions, where x and k± are the longitudinal momentum fraction and 



11 



transverse momentum of the quark, respectively. For a longitudinally-polarized nucleon, 
one introduces a longitudinally-polarized quark distribution gn(x,k±) = Aqi(x,kj_) and 
a transversely-polarized distribution hj~ T (x,k±) = 5qi{x, k±). Finally, for a transversely- 
polarized nucleon, one introduces a spin-independent distribution f^ r (x,k±) = qT(x,k±) 
arising from final and/or initial state interactions, and a longitudinally-polarized quark 
polarization gir{x,kj_) = Aq T (x, k±), a symmetrical transversely-polarized quark distribu- 
tion hiTiXyk^) = SqT{x,k±) and an asymmetric transversely-polarized quark distribution 
hj~ T (x, k±) = 5qT'(x, k±). 

Out of the eight distributions, four of them are related to the nucleon helicity flip and are 
the main interest in this subsection: unpolarized and longitudinally-polarized quark distribu- 
tions in a transversely-polarized target, qr{x, k±) and AqT(x, k±), and transversely-polarized 
quark distributions in an unpolarized and longitudinally-polarized nucleon, 5q(x, kj_) and 
5qi(x, k±). Two of them, qr(x,k±) and Sq(x, k±), are non-zero only when there are initial 
and final state interactions. 

Unlike the usual Feynman parton distributions, there are some important subtleties about 
the transverse-momentum dependent parton distributions. First is gauge invariance and uni- 
versality. Each of these distributions can be gauge-invariantized in different ways depending 
on choices of gauge links and the results are different. Which one appears in a particular 
process requires a careful study. In the Drell-Yan process studied by Collins and Soper ||8| , 



the distributions are defined in axial gauges. In deep-inelastic scattering, it has been shown 
that the gauge links are in the future direction along the light-cone. If the light-cone gauge 
is used, an additional gauge link at spatial infinity is required to render the distributions 
independent of additional gauge fixing fi3]| . Different versions of the distributions appear 
in different hard processes imply that the universality of these distributions are lost. For 
example, the transverse-momentum dependent parton distributions measured in DIS cannot 
be used to make predictions of the cross sections for Drell-Yan. For processes such as jet 
production in hadron-hadron scattering, it is not yet clear which version of the distributions 
is relevant. 

Another issue is infrared divergences. The distributions involving light-cone gauge links 
have infrared divergences, which are cancelled after integrating out parton transverse mo- 
mentum. Physically, the infrared divergences can be understood in the following way. The 
transverse-momentum dependent distributions can be measured in DIS through single-jet 
production. However, the single-jet cross section is not well-defined because the vertex cor- 
rection has infrared divergence. The divergence is cancelled by the soft gluon radiation. 
Since a hard gluon radiation is considered as a two-jet event, one must introduce an infrared 
cut-off to separate the single from double jet event. In the following discussion, we implicitly 
assume the infrared cut-off is there when needed. 

We now use the wave function of the previous section to calculate the helicity-flip parton 
distributions. We first consider the case of a transversely polarized target. Introduce the 
quark density matrix, 

M = p + J ^^e^'-^\PS\mLlL^(0\PS) , (43) 

where is the gauge link along the light-cone in the covariant gauges, S M is the polarization 
vector of the nucleon normalized to S^S^ = — 1, is a light-cone vector such that p~ = 0. 
The distributions we are interested can be obtained from M. through the expansion ||22|| , 

M = m \ qr{x ' k ^ tilVa ^^uk a S fi + Aq T (x, k ± ) l5 j(k ± ■ Sj_) + ...] , (44) 
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where M is the nucleoli mass. Inverting the above equation, we obtain 

&q T (x,k ± ) = JL-J ^||ie^ + ^-^^)(P^|^(0)Lj 7 + 75 L^(OI^±) • (45) 
where the transversely-polarized nucleon in the x direction is related to the helicity states 

by \S X ) = (| T) + | i))/V2. 

If these distributions are calculated in the light-cone gauge with the advanced boundary 
condition for the gauge potential, the gauge links can be ignored, but the wave function 
amplitudes are complex jEJ. Express the quark fields in term of Fock operators (x > 0), 

A 

where V3 is the 3-space volume. Using the above and the proton wave function from the 
previous section, we obtain 

Aq T (x,k±) = J d[l]d[2]d[3]Vx 1 2x 2 2x 3 Re[F q ] , 

qr(x,k ± ) = J d[l]d[2]d[3Wx 1 2x 2 2x 3 Im[F q ] . (47) 

The functions F q for the u-quark is 

F u = 2 [8®(k - fc0^ 2) *(l,2,3)^f 4) (1,2,3) - 6®(k - h)^ 2 \l, 2, 3)$ M) *(2, 1,3) 

+{S ®(k - h) + 6®(k - fc 3 ))^ 1>2) *(l, 2, 3)(^ (3 ' 4) (2, 1, 3) + ^ (3 ' 4) (2, 3, 1)) 

+(6®{k - h) + 5 (3) (fc - k 2 ))^->(l, 2, 3)^ 6+ )(l, 2, 3) } , (48) 

where 

<y(3)(fc _ ki ) = 5 (x - x t )5 (2 \k ± - k i± ) , 
$ 1,2) (1,2,3) = ^(1,2,3)^ -(^-A;?^«(l,2 ) 3))A ; i e«, 
^(1,2,3) = M^(l,2,3) + ^#)(1,2,3) , 
^ 5 -)(l,2,3) = ^(1,2,3) -^(1,3,2)) , 

^ 6+ )(l,2,3) = -k 2 l± k 3L • k ± ^ (1,2,3) + (k 2 J 3± ■ k ± + k 2 3 J 2± ■ k ± )^ 6 \2, 1,3) 

+k 2 1± k 2 -k ± (^\l,3,2) + ^ (6) (2,3,1)) . (49) 

For the d-quark, on the other hand, we obtain 

-fc 3 )^ (1 ' 2) *(l,2,3)^ {3 ' 4) (l,2,3) 

-5® (k - kSf' 2 \l, 2, 3)(^ (3 ' 4) *(2, 1, 3) + ^ (3 ' 4) *(2, 3, 1)) 
+5®(A;-A;3)^ 5 - ) *(1,2,3)^ 6+ )(1,2,3)} . (50) 
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If all the wave function amplitudes are real, qT(x,k±_) vanishes identically. qT(x,k±_) was 
first introduce by Sivers [Jl9| , |21|j and it characterizes the azimuthal asymmetry in the quark 
transverse momentum distribution when the nucleon is transversely polarized. There is 
much interest in the this distribution in interpreting the single spin asymmetry measured in 
electron-proton deep inelastic scattering ^ |5D], [H], [53], . 

The other two helicity-flip distributions can be obtained from the quark density matrix 
M through the following projection 



M = 2M k± ^ k ^ + 6qL ( x > k ±) ia ^^pMS ■ n) + ...] . (51) 

The transversely-polarized quark distribution in an unpolarized nucleon Sq(x, k±) is novel, 
and has the similar physical origin as that of the experimental phenomenon where a hy- 
peron produced in unpolarized elastic proton-proton scattering is polarized in the direction 
perpendicular to the production plane [[55 , |56[| . Invert the above equation, we arrive at 

Sq L (x,k ± ) = ^yj ^|^e i ( fe+r -^^ ) (P5|#(0)4 ( 7 + ^(e)|P^|) • (52) 

where % can take either x or y. 

Inserting the plane wave expansion for the quark fields, we get 

W*>k ± ) = 4x{ l Yv j (M( k ^( k ) -q}(kMk)\P) , 
SQL(x,k ± ) = ™ k . M ^ )iv J (PSM(k)9r(k) - q\(k) qi (k)\PS\\) ■ (53) 
Using the above and the proton wave function from the previous section, we obtain, 
Sq L (x,k ± ) = y d[l]d[2\d[3Wx 1 2x 2 2x 3 Re{H q ] , 

5q{x,k ± ) = rj f d[l]d[2]d[3}Vx 1 2x 2 2x 3 Im{H g } . (54) 
k±_ J 

The function H q is 



H u = 2 {-^(A; - ^^'^'(l, 2, 3)(^f' 4) *(3, 1, 2) + ^ 3 ' 4) *(3, 2, 1)) 

-5®(k - k 2 )(i>f' 2) ' (1, 3, 2) + ^ (1 ' 2) '(2, 3, 1))^ M) *(1, 2, 3) 

-5^(k - 2, 3)(^( 5 ~)(2, 1, 3) + k{ft 5 -\l, 2, 3)) 

+5®{k - h)(ft 6+ >(l, 2, 3)^(1, 2, 3) + ft 6+ >(2, 1, 3)^(1, 2, 3)) } , (55) 

where $ 1,2) '(1, 2, 3) = ^ (1) (1, 2, 3)k j + (kfk^ - k\kf)^{l, 2, 3))A;V 3 '. For the down quark, 
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H d = 2{5^(k-k 3 )^ (1,2,3)^(1,2,3) 

-S^(k - k 2 )^' 2) *{l, 2, 3)(k{^ 5 -\l, 2, 3) + fc^( 5 ")(3, 2, 1)) 

+<jW(* - fcO^+'^l, 2, 3)^(1, 2, 3) - ^ 6+ )*(2, 3, 1)^(1, 2, 3)) } , (56) 
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where 

^ 6+ '>(l,2,3) = (fc 2 \4±-fc± + ^L^±-^)(^ (6) (2,l,3) + ^ (6) (3,l,2)) 

+k\ ± k 2 ■ fc^ (6) (2, 3, 1) + k 2 1± k 3 ■ fc ± ^ (6) (3, 2, 1) . (57) 

If all amplitudes are real, Sq(x, k±) vanishes identically. 
B. Twist-Three Parton Distributions 

It has been known for many years that the helicity flip phenomena in inclusive scattering 
are of higher twist effects (twist-three, to be exact). For example, the single spin asymmetry 
measured in pion production in polarized baryon (proton, hyperons) and unpolarized target 
scattering arises from hadron helicity flip J57|. The asymmetry is expected to vanish asymp- 
totically like 1/kx, where k± is the pion transverse momentum ||58|| . The other example 



is the (?2 structure function measured in inclusive DIS with transversely polarized nucleon 
targets. The structure function arises from the interference between photon scatterings with 
longitudinal and transverse polarizations and involves the helicity flip of the nucleon target 
because of the angular momentum conservation p5| . The associated asymmetry goes like 
1/Q where Q is the virtual photon mass. 

In both examples, the baryon helicity flip must be reflected in the hard scattering sub- 
processes. However, if the quark transverse momentum is neglected as one usually does in 
the Feynman parton model, the helicity flip cannot be managed through a massless quark; 
this has been considered as one of the difficulties in perturbative QCD (|Tj^, It turns 
out that the leading contribution to the baryon helicity flip comes from the quark orbital 
angular momentum and associated transversely-polarized gluons El]]. This can be seen 
from the following light-cone expression for the g 2 (x) structure function 



9t(x) = ^e ax (PS ± \m fl l5 $ ± i P ± (\n)^(Xn)\PS ± ) + h.c. . (58) 

In our approximation, the gluon potential in the covariant derivative is neglected. Using the 
three-quark wave function, we find, 



9t(x) = J kiAq T (x,k ± )d 2 k ± . (59) 

In fact, the above result holds for any nucleon wave function so long as the gluon potential 
in the covariant derivative can be ignored ||23|| . A word of caution, however, is appropriate. 
When the gauge potential is neglected, gr(x) is no longer gauge invariant. How it is possible 
then to express gr{x) in terms of the the gauge-invariant Ag^(x, k±)7 The answer is that 
the above relation is only true in light-cone gauge with the advanced prescription for the 
light-cone singularity. In any other gauge, the two have no relation. 

It is interesting to note that the rotational invariance imposes the following condition on 

9t(x) m 

.1 ,1 

dxgr(x) = / dxg±(x) . (60) 
'o Jo 

Therefore, when the gluons are neglected, the rotational symmetry demands the quarks have 
non-vanishing orbital angular momentum. This is indeed true for a massless quark in the 
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MIT bag model. However, it is possible that gr(x) may have a delta function at x = [62 
In this case, the above is useless in constraining gr(x) at any non-zero x. 

Similar analysis can be carried out for the twist-three distribution hi(x) which has the 
following expression, 

h ^ = T~m I ^e lXx (PS ± \m P ± (\n)^(Xn)\PS ± ) + h.c. . (61) 

4xM J Z7T 

Using the proton wave function, we find 

h L {x) = - j k 2 ± 6q L (x,k L )d 2 k ± , (62) 

which is again consistent with [23] apart from a factor of 2. The twist-three distribution 
Hl(x) can be measured in Drell-Yan scattering with the transversely-polarized protons scat- 
tering on longitudinally-polarized protons pi. 



C. Generalized Parton Distributions 

Generalized parton distributions (GPDs) were originally introduced to describe the an- 
gular momentum structure of the nucleon [Q, |B5| . D. Mueller et al. encountered the same 
distributions earlier in searching for an observable whose scale evolution interpolates be- 
tween those of the Feynman parton distributions and the quark distribution amplitudes in 
mesons [f|]. Radyushkin defined the so-called double distributions which contain essentially 
the same information The new distributions contain all the form factors of the twist- 
two operators || and, in a special kinematic limit, reduce to the Feynman distributions. In 
this subsection, we calculate the helicity-flip GPDs using the three-quark wave function of 
the proton. 

At the leading twist, there are eight GPDs for each quark flavor |)5| |66|] , out of which four 
involve the nucleon helicity flip. In this paper we consider only E(x, £, Q 2 ) which is directly 
relevant for the spin structure of the nucleon fl6j . Extending the calculation to other three 
is straightforward and will not be presented here. The definition of E(x, £, Q 2 ) follows from 



dX , 



2tt 



P' 



1 1 in^n A 

H q {x^Q 2 ) -U(P') iU{P)+E q {x^Q 2 ) -U(P') 2 £ » U{P) 



(63) 



where the explanation of various notations can be found in the original references. 

To evaluate E(x, £, Q 2 ), we choose a coordinate system in which the momentum transfer 
vanishes along the + direction, 



0, Q 2 



(64) 



Then the initial and final nucleons have the same P + momenta. Their transverse momenta 
are equal in magnitude and opposite in direction (a Breit frame in the transverse direction). 
The energy- momentum conservation constrains q~ = 0. With these choices, we can only 
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calculate the distribution at the skewness variable £ = 0. [For £ 7^ 0, one can use the 
coordinates in Ref. |j7| .] 

Our calculation follows a similar formula for F 2 (Q 2 ) in Ref. [ |TC|j . 



(65) 



For x > 0, the bilocal operator reduces to quark creation and annihilation operators in the 
momentum space, and for x < 0, it becomes antiquark creation and annihilation. The later 
contribution vanishes in the valence approximation. 

Using the proton wave function from the previous section, one finds 

£ M (x,£ = 0,Q 2 ) = — — - / d[l]d[2]d[3] V*i2x 2 2x 3 {A^(/>^(1, 2, 3) 
— q x + iqy J 

-A (3 ' 4) (1,2) '(l,2,3) -A (5) (6) (l,2,3) + A (6) (5) (l>2,3) } , (66) 

where 

</> (1 ' 2) (l,2,3) = ^(l,2,3) + z(^|-A;M^ 2 )(l,2,3), 
^'(1,2,3) = ^ 1 )(1,2,3)-^|-A;^)^ 2 )(1,2,3), 
(3 ' 4) (1,2,3) = A;r^ 3 )(l,2,3) + A; 2 -^ 4 )(1,2,3) , 
(5) (1,2,3) = fc+(^(l,2,3)- ^(1,3,2)) , 

(6) (1,2,3) = k^k^ 6 \l,2,3) - k^k^ (1,3,2) . (67) 
And the functions A are defined through 

A (1,2) = 8{x - xx)0 (1 ' 2) *(2', 1", 3') + S(x - x 2 ){2^ 2 >{2", l', 3') + cj>^ 2 >(3', l', 2")) 
+6{x - x 3 )(0 (1 ' 2) *(2', 1', 3") + V^ 2 )*(3", 1', 2')) , 
= 5{x - x 2 )0 (3 ' 4) *(2", 1', 3') + 6(x - x 1 )(20( 3 ' 4 )*(2', l", 3') + <^*(2', 3', l")) 

+5(x - x 3 )(</) (3 ' 4) *(2', 1', 3") + ^ (3 ' 4) *(2', 3", 1')) , 
= 5(x-xi)(0 (5) *(l",2',3') + (5) *(2 / ,l /, ,3')) 
+5(x - x 2 ){^>{\\ 2", 3') + 0< 5 )*(2", 1', 3')) , 
A® = ^(x-s 1 )(0 (6) *(l // ) 2 / ,3 / ) + (6) *(2^1^3 , )) 

+5(x - £ 2 )(0 (6) *(1', 2", 3') + <P^*(2\ 1', 3')) , (68) 

where the transverse coordinates are %' = k{ — Xiq± and i" — k{ + (1 — Xi)q±, kf = kf ± i/cf . 
The result depends on the interference of the wave function amplitudes with different quark 
orbital angular momentum. 
Similarly for the d-quark, 

o/i/f r 

E d (x^ = 0,Q 2 ) = ——— / d[l]d[2]d[3Wx 1 2x 2 2x 3 {B^ ( j ) ^\l,2,3) 
—q x + iqy J 

_£(3,4)^(l,2)' (1; 2; 3) _ B (5)^(6) (1; 2; 3) + 5 (6)0(5) (1; 2 , 3) } , (69) 

where 

= 5(x - x 3 )0 (1 ' 2) *(2', 1', 3") + S(x - x 1 )(0( 1 - 2 )*(2 / , 1", 3') + <^ 2 >(3', 1", 2')) , 
= 5(x - x 3 )0 (3 ' 4) *(2', 1', 3") + 5(x - x 2 )(0 (3 ' 4) *(2", 1', 3') + 0( 3 ' 4 >*(2", 3', 1')) , 
£(5) = 5(x-x 3 )(0 (5) *(l , ,2 , ,3") + (5) *(2 / ,l / ,3")) , 

B(6) = 5( x - x 3 )(0 (6) *(l', 2', 3") + (6) *(2', 1', 3")) • (70) 
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Although we have allowed q to be arbitrary in the above formulas, it is simpler to take it in 
the x direction. With this choice, it is not difficult to see that the distributions are real. 



D. The Pauli Form Factor 

The F 2 form factor can be obtained from the sum rule ||, 

F2(Q 2 ) = j dx[e u E u (x, C Q 2 ) + e d E d (x, £, Q 2 )} , (71) 

where we have neglected the strange quark contribution. Recently, the Pauli form factor of 
the proton has been extracted from the recoil polarization [§, ||, which differs significantly 
from the previous extraction. In the valence approximation, independent of Q 2 , a non-zero 
^{Q 2 ) is an indication of non-zero orbital angular momentum of the quarks in light-cone 
quantization. This point has been stressed in a number of recent papers which seek to 
explain the new Jlab data pR 



In large-Q 2 limit, we can use the asymptotic behavior of amplitudes to derive the asymp- 
totic behavior of F 2 (Q 2 ). The result will be related to the leading and twist-three light-cone 
amplitudes introduced by Braun et al. @]. We leave this subject to a future publication. 

Finally, we pointed out that the S-form factor of the gravitational current can be ob- 
tained from the second moment of E(x,£,,Q 2 ) [[| and is also helicity-flipping. In Ref. |I2| , 



the B form factor is calculated for the electron in quantum electrodynamics, with explicit 
involvement of orbital angular momentum. 



V. SUMMARY AND CONCLUSION 

This paper is motivated by determining the components of the proton wave function 
with non-zero orbital angular momentum. This of course can only be down in a certain 
truncation scheme, because the field-theoretical nucleon has an infinite number of non- 
vanishing amplitudes and there is no finite set of experiments which can determine them 
completely. 

We use the valence quark models as a guide to truncate the light-cone expansion of the 
proton wave function. We can view this either as a starting point for a systematic expansion 
or as an effective theory in which gluons and sea quarks have been "integrated out" . 

We classify the light-cone amplitudes with three quarks, and find that six independent 
amplitudes are needed for a complete description: two with L z = 0, three with L z = 1, 
and one with L z = 2. If the light-cone gauge fixing is invariant under combined parity and 
time-reversal, these amplitudes are real. Otherwise, they are complex. Apart from one of 
the L z = amplitudes, the other five contain nontrivial quark orbital motion. 

We calculated a number of nuclear observables which depend at least linearly on the 
amplitudes with non-zero orbital angular momentum. Therefore, a non-vanishing result of 
the experimental measurement is an unambiguous indication of the quark orbital angular 
momentum, barring the neglected gluonic contributions mentioned above. 

It is our hope that multiple observables of this type can be explained by a unified set 
of phenomenological amplitudes with nontrivial orbital angular momentum. If not, one 
may systematically go beyond the minimal Fock component by including, for example, the 
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gluon contributions. Ultimately, we hope to have a semi-realistic picture about the intrinsic 
deformation of the proton from experimental data. 

The authors thank S. Brodsky for a number of useful discussions. This work was sup- 
ported by the U. S. Department of Energy via grants DE-FG02-93ER-40762. 
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